Nonlinearity effects in the kicked oscillator 
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The quantum kicked oscillator is known to display a re- 
markable richness of dynamical behaviour, from ballistic 
spreading to dynamical localization. Here we investigate the 
effects of a Gross Pitaevskii nonlinearity on quantum mo- 
tion, and provide evidence that the qualitative features de- 
pend strongly on the parameters of the system. 

PACS numbers:05.45.-a 



The dynamical behavior of quantized area-preserving 
maps has proved to be one the most relevant field in the 
discipline of quantum chaos (see Q): in particular the 
discovery of quantum dynamical localization j^], namely 
the quantal suppression of classical deterministic diffu- 
sion, has provoked a vast amount of theoretical and ex- 
perimental work. The paradigmatic example in this field 
is the quantum kicked rotator (see ||), obtained upon 
quantization of the classical standard map [Q. From a 
classical point of view the system is of KAM nature: for 
small values of the stochasticity parameter global trans- 
port is inhibited by invariant curves: once the last in- 
variant torus is destroyed, transport properties abruptly 
change, and, in typical situations of strong chaos, there is 
a diffusive spreading in momentum |IJ|5|] , which character- 
izes also the quantum motion for times shorter than the 
break-time if,, where quantum localization regime sets in 
and the momentum spreading is suppressed Jg]. We re- 
mark that this picture is valid for generic values of the 
effective Planck's constant: quantum resonant motion, 
characterized by ballistic spreading, appears when h as- 
sumes rational values 0,0 • 

Another example of a quantum system originating 
from a two dimensional area preserving map is the kicked 
harmonic oscillator (see |p|-pl]| and references therein): 
the classical hamiltonian is 

H{p,x,t) = — i-p 2 + ^oj 2 x 2 + ecos(k x)5 To {t) (1) 

where the time dependence is through the periodic delta 
function 
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In particular the resonant case (T = 2itp/q) is charac- 
terized by the presence of a stochastic web (for arbitrar- 
ily small values of /C) supporting unbounded transport 
]l2| , p~3|] , while in the non-resonant case a threshold JC{E ) 
exists below which unbounded transport is not sustained 
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The kicked harmonic oscillator has been proposed as a 
model of different physical phenomena: from electronic 
transport in semiconductor superlattices |l5f| , to ion traps 
Pq |. In the latter case the harmonic potential is repre- 
sentative of the ion trap, while the kicking term arises 
from a time periodic standing wave laser field. Obvi- 
ously such examples require a proper quantum mechani- 
cal treatment of the Hamiltonian ([l]), the corresponding 
Schroedinger equation being (once expressed in dimen- 
sionless variables t = u>ot and x = \/ mouo/hx) 



2dx 2 



so that quantum dynamics depends upon three parame- 
ters 
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T = T q lo q (4) 



By rescaling variables x = kgx and t — toot we realize 
that dynamics is only dependent on the parameters 



Once again the behavior is quite sensitive to number the- 
oretic properties of T: in particular the crystal cases (h| 
T = 2n/q with q S {1,2,3,4,6} admit a one-parameter 
group of commuting generalized translations commuting 
with the hamiltonian (exceptional parameter values JlCj ] 
may also lead to two-parameter groups): the correspond- 
ing dynamical behavior is diffusive (or ballistic in the 
exceptional cases). We remark that the q = 4 case cor- 
responds to the (symmetric) kicked Harper model [ p"7| |. 
Outside resonant parameter values there are indication of 
a localization-delocalization transition for resonant-non 
crystal cases tyj, while the simulations reported in |1C] 
for non-resonant cases suggest dynamical localization | Pc ] 
(however we observed a derealization transition in the 
irrational case too). 

Recently it has been suggested [|l9| how the widespread 
interest and experimental activity in Bose-Einstein con- 
densation |2(J makes it natural to study the effect on 
Gross-Pitaevskii nonlinearities on the kicked oscil- 
lator (thus turning it into a model of a trapped con- 
densate under a laser field in the spirit of [|l6|). The 
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Gross-Pitaevskii nonlinear correction to the Schroedinger 
equation is of the form u|7/j| 2 ?/>, where the coefficient u is 
of the same sign of the scattering length p3] (we will 
mainly deal with a positive u in what follows): using the 
same rescaling in dimensionless variables mentioned in 
the quantum case the equation reads 
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Even if the cubic nonlinearity acts like an effective repul- 
sive potential, the main observation in as regards 
the dynamical effect of the Gross-Pitaevskii nonlinearity 
in a crystal q = 6 case, has been its tendency to oppose 
quantum spreading: it was suggested that this is due to 
a breakup of quantum symmetries for nonzero nonlinear- 
ity. Before presenting the results of our simulations we 
have to mention that nonlinearity effects have also been 
considered a few years ago for the kicked rotator [ p3|j24| 
(for a cubic nonlinearity of opposite sign). Here the sce- 
nario is quite different: when the nonlinearity is absent 
the system exhibits quantum dynamical localization: a 
sufficiently strong nonlinearity may then induce chaotic 
transitions between localized modes, leading to (subdif- 
fusive) derealization. 
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FIG. 1. Classical transport along the stochastic web. 

To investigate the effect of the nonlinearity, we stud- 
ied the equation (|5|) in two different regimes: a crystal 
(q = 4) case, and an irrational case (T = tt/(t/5 + I)). 
The evolution of the the kicked oscillator is conveniently 
studied by using a discretized propagator Eg]: 



ij){x',t') = / dxG(x',x;t)ip(x,0) 



(7) 



where 



G{x',x;t) = Cexpl ' m ° / a, ° N ((x 2 + x' 2 ) cos(^) - 
l_ 2n,sm(ti;oi) 

(8) 

and once the discretized positions Xi = (i — N/2)A X are 
introduced we have that the propagator Q is unitary if 
we put 
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In the coordinate representation the action of kicks is 
multiplicative on the wave function. 
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FIG. 2. Quantum transport (g = 4 one parameter symme- 
try group). See the text for details. 

Simulation of quantum evolution is considerably more 
complicated once we introduce a nonzero nonlinearity: to 
propagate the wave function between kicks we separate 
the time independent part of the hamiltonian into the os- 
cillator and the nonlinear part, and use the lowest order 
split method p5|] (this typically requires using ten time 
steps between consecutive kicks in order to get stable re- 
sults). As initial state we consider the ground state of 
the GP equation (without kicks) shifted in in the chaotic 
region nearest to the origin. The ground states for differ- 
ent values of the nonlinearity parameter arc obtained by 
evolving an eigenstate of the quantum harmonic oscilla- 
tor under the imaginary time Gross-Pitaevskii equation 
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The first case we take into account is a crystal q = 4 
example. In fig. (1) we show the classical phase pic- 
ture, exhibiting unbounded transport along the stochas- 
tic web. We put T = ir/2 (and ujq fixed in such a way that 
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generalized translations form a one parameter group of 
symmetries), e = 0.7 (all others linear quantum parame- 
ters are fixed by taking ft = 1 and £ = y/2: we will adopt 
this choice for other examples too). 

This case is characterized by a one-parameter group 
of simmetries (generalized phase space translations) and 
thus the quantum case is expected [ jlCfl to exhibit a dif- 
fusive momentum spreading (the evolution corresponds 
of the upper line in fig. (2)): the effect of nonlinearity 
is considering by taking u — 0.5, 1, 2.5 and 5: the corre- 
sponding curves are shown in fig. (2). Such simulations 
are performed by using a N = 2 14 discretization of the 
position variable (which has been checked to provide a 
reliable choice up to the considered evolution time, by 
comparing the results with a simulation with twice the 
number of points). To smooth out oscillations in the 
evolving patterns we plot the integrated second moment 
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FIG. 4. Quantum transport (q — 4 two-parameter symme- 
try group) group, e = 0.7. The configuration space has been 
discretized with 2 12 points. 
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FIG. 3. P 2 (t), for u = +0.5 and u = -0.5, and same 
conditions as fig. (2). 

The qualitative features confirm the observation in 
Jl]| ) namely that the most striking effect of the non- 
linearity is to oppose quantum delocalization. This has 
been claimed to be due to symmetry breaking effects of 
the nonlinearity, inhibiting transport along delocalized 
Floquet states. We remark that also at a classical level 
related features have been observed: if noise is added to 
the kicked Harper map, transport along the stochastic 
web is slowed down p7| . 



In principle a positive nonlinearity acts like a repul- 
sive potential, but here the symmetry breaking effect is 
not related to the sign of the effective potential, as we 
see from fig. (3), where it is shown how the sign of the 
nonlinearity has a tiny effect on momentum spreading. 
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FIG. 5. Quantum transport irrational case. Discretization 
is over 2 14 points. The inset shows the kicked oscillator case 
on a longer time scale. 

If such a view is correct, the effect must equally ap- 
pear in the resonant crystal case, where a group of two- 
parameter symmetries leads to ballistic transport flio| |: 
such an effect is indeed evident even for short times, see 
fig. (4): after a characteristic time scale, which shrinks 
as nonlinearity increases, the deviation from the kicked 
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oscillator case is more and more marked as u increases. 

A priori the situation is different when we consider the 
oscillator outside the crystal regime: to this end we ana- 
lyzed the case in which e = 1.4 and T = n/(y/5 + 1). In 
this case the kicked oscillator displays dynamical local- 
ization: the striking observation is that here the nonlin- 
earity acts in an opposite way, enhancing the quantum 
derealization, see fig. (5). So, when symmetries are not 
present in the quantum case, nonlinearity seems to play 
a completely different role. This is at least in qualitative 
agreement to what happens to the kicked rotator evolving 
under a nonlinear Schroedinger equation, or even when 
noise is superimposed to the quantum evolution 28 1 . We 
have checked that the same happens even for higher val- 
ues of e, when the oscillator undergoes a derealization 
transition. 

In conclusion, we have analyzed how nonlinearity in- 
fluences a complex, and physically relevant quantum sys- 
tem, the kicked harmonic oscillator. We provide evi- 
dences that, at least at moderate times, it opposes quan- 
tum diffusion when transport is linked to symmetry prop- 
erties of the linear hamiltonian, while it may lead to dif- 
fusion enhancement when no symmetry breaking occurs. 
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